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ABSTRACT 


This study investigates the maneuvering of ships in- 
volved in the replenishment at sea operation under calm 
water conditions. 

Two sets of linear differential equations of motion 
of a ship in three degrees of freedom are implemented for 
an analog-digital simulation (hybrid operation). Mainly 
a two phase hybrid simulation is carried out. 

In the first part the real-time dynamic response of 
a Single ship is obtained. Small perturbations, result- 
ing from small changes in rudder ae and propeller RPM, 
are studied. 

In the second part interaction forces and moments 
are applied only to the leading ship as it is overtaken 
by the tracking ship. Thus the response of the leading 
ship is of primary interest on this phase. 

Finally interaction forces and moments are applied 
to both leading and tracking ship. Thus the complete 


response of the UNREP operation is obtained. 





EL. 


ET. 


IV. 


TABLE OF CONTENTS 


Нино Е == ССС А 
INTERACTION EFFECTS- -— - - - - - - - - – – — 
DERIVATION OF THE LINEAR MATHEMATICAL MODEL- 
A. GENERAL CASE-SIX DEGREES OF FREEDOM- 一 一 
B. HORIZONTAL PLANE MOTION— — — — — — — — — 


C. LINEARIZATION THROUGH TAYLOR'S 
SERIES EXPANSION — — — — — — — — — — — — 


ПИОН ЫШ МЕ КЫ О PPEQUATION = = — — — — — — 


AXES FIXED RELATIVE TO EARTH—LANGRANGIAN 
SYSTEM — — — — — — — — — — — — — — — — — — — 


ANALOG COMPUTER PROGRAMMING FOR EACH 
SHIP INDEPENDENTLY ~ = = = = — — — — — — — — 


A. PRELIMINARY STUDY— — — — — — — — — — — — 
1. Linearized Equations of Motion — — — 


2. Magnitude Scaling of the 
Pincaelzedwbguaslens = - — = — 一 一 一 


3. Calculation of Potentiometer 
Coefficients ------------ 


I "naNesEParehdggeeonfiguratron- = = = 
a. Static Test— — — — — — — — — — 一 


B. CHARACTERISTIC LINEAR RESPONSE OF 
MARINER= = — — — — — — — — — — — — — — — 


l. Dynamic Test — — — — — — — — — — — — 
2. Steering Control — — — — — — — — — — 
AA ys ls Of Steering Control 一 一 


b. Time Lag and Settling Time of 
Imoroved Steering Control= = = — 


10 
13 
13 


15 


16 


19 


24 


25 
25 


25 


26 


28 
36 


36 


41 
45 
33 


55 


60 





WI, 


RN LON OF COORDINATES AXES—- = = = 61 


l. Introductory Discussion ------- 61 

2. Computed Linear Response of Mariner 
after Coordinates Transformation- — — 76 
INTERACTION FORCES AND MOMENTS ARE INCLUDED — 92 
A. ΥΕ DISCUSSION — — — — — — — = — 92 


BEZZINEERAECTTONZBATA CURVES— — 
MERCO CATTON PROGRAMMING = = = =-=- — — — DES 


C. MEASUREMENT OF LONGITUDINAL AND LATERAL 


ВИКИ ОГ шотаАСЕк5——= = = =— = ーー ーー — — 96 
l. Longitudinal Distance — — — — — — — — 96 
2. Lateral Distance- = =- =- =- == == -= - 97 


D. APPLICATION OF INTERACTION MOMENT N AND 


FORCE Y TO SHIP A (LEADING SHIP)- -— — — — 100 
ШОШО оеше зову DISCUSSJONH = = = === = 100 
2. Dynamical Representation of Inter— 
action Moment and Force — Analog 
A M το ο - ТО? 
3. Obtained Responses for Phase I- — — — 105 


a. Stationary Runs (i.e. both ships 
have same propeller speed)- — — — 105 


b. .Runs with Different Propeller 
Speeds between Ship A (Leading 
Snip and Ship B (Tracking Ship)— 126 


4. Obtained Responses for Phase II - — — 135 


a. Stationary Runs (i.e. both ships 
have same propeller speed)— — — — 148 


b. Runs with Different Propeller 
Speeds between Ship A (Leading 
Sito ана өлір В (Туаскіпа Ship)- 148 


С. Manual Control is Applied in 
Phase II— — — — — — — — — — — — — 175 





PO ου το TONZZZREMARKS = SUMMARY — =- = = -= 


INVESTIGATION OF THE CONTROL PROBLEM — — — — 


S I TS DIII TIC PTOODP= —"— — — =— = = — — — 


PERDIS STAN ERORTSTATION KEEPING LOOP — — — — 


COMPUTER PROGRAM — — = — — — — — — — — — — — — — — 


LIST OF REFERENCES — — — — — — — — — — — — — — — — 


E TEE ION Он = Е а τ ΞΞΞΞΕΞ 一 


FORM DD 1473 


4a 


188 


196 


235 


236 


231 





р,а,г 


e 


TABLE OF SYMBOLS 


SP ojj Criterion for stability in 
Straight line motion 


Moments of inertia about the X, Y, Z axes, 
respectively 


Rolling, pitching, and yawing moments, 
respectively 


Ship length between perpendiculars (LBP) 
Mass of ship 


Angular velocities of roll, pitch, and yaw, 
respectively 


Velocity of the origin of the body axes 
relative to the fluid 


Time and time interval, respectively 


Velocity components of the origin of the 
body axes relative to the fluid (longi- 
tudinal, and normal components, respective- 


ly 


Initial equilibrium velocity component 
(ahead straight line motion at constant 
speed with rudder at amidships) 


Uu — Us 

Acceleration components of the origin of 
the body axes relative to the fluid (longi- 
tudinal, transverse, and normal components, 
respectively 


Hydrodynamic force components in ship body 
(longitudinal, lateral and normal components) 


Coordinate axes fixed in ship. Origin of 
axes system need not be at the center of 
gravity of the ship (positive direction 
forward, starboard, and downward, respec— 
tively 





ф,9,% 


An 


AR 


Geoeccainates ot the center of mass of the 
ship relative to body axes 


Coordinates relative to the coordinate 
system fixed in the earth 


Coordinates of the center of mass of the 
ship relative to the coordinate system 
fixed in the earth 


Zuelesor drift. Lateral side to side 
Separation distance between ships 


Angular displacement of a control surface, 
normal to the rudder angle 


Angles of roll, pitch, and yaw, respec— 
tively 


Mass density (mass density if sea water is 
1.9905 lb-sec 


Velocity components of the ship axes 
relative to the space coordinate system 


Longitudinal separation distance between 
midships 


Small change in RPM propeller speed 


Small change in rudder angle 





ACKNOWLEDGMENT 


The author wishes to express his sincere appreciation 
to Dr. George J. Thaler for guidance, assistance and con- 
tinuous encouragement which he provided during the pursuit 
Bernie study. The author would also like to express his 
appreciation to Reidar Alvestad, Samuel H. Brown, Bob 
Limes and Lieutenant David Hall for their valuable assist- 
ance and information given to him. Finally the author 
would wish to thank his wife Niki for her assistance in 
typing this study as well as in drawing the various figures. 
The author is also deeply indebted to A. G., the great man 


who made it all possible. 





ie tT RODUCTION 


At the present time the operational procedure of re- 
plenishment at sea underway is described by two manuals 
of the Department of the Navy [1] and [2]. This procedure 
includes collision hazard. Therefore an analysis of the 
replenishment operations may provide a mean for improve- 
ment of ship maneuvering and ship control. The maneuver 
of replenishment at sea involves six factors: course, 
speed, distance between ships, the approach, station keep- 
ing and departure. 

This thesis describes the E of underway replen- 
ishment (UNREP) using an analog computer simulation. This 
analog computer simulation model will describe two Mariner 
Class merchant ships. An effort is made to define the 
control parameters necessary to provide information for 
the conning officer and/or helmsman which will increase 
the safety conditions of the UNREP operation. According 
to the tactical requirements the UNREP operation is inter— 
preted as follows: 


a. The replenishing (leading) ship is responsible 
for course keeping only. 


b. The receiving (tracking) ship is responsible for 
both course and station (distance between ships) 
keeping. 


For this reason one ship is constrained to move in a 


straight line course at constant speed. The rudder and 





propeller of the second ship are controlled in order to 
keep it in a parallel course relative to the first ship. 
The analog computer output is presented and analyzed 
to determine the necessary information for the ship con- 
trol parameters required during these simulated UNREP 


operations. 





INE ACTION TEE FECTS 


oen underway, the pressure distribution on the ship 
hull varies (as is shown in Figure 1) due to the variable 


velocity of the water flow, the so-called venturi effect. 
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Figure 1. Pressure distribution on the hull. 


When two ships are alongside underway, this venturi 
effect is increased and becomes more complicated due to 
the mixing of the pressure areas of two ships. 

Figure 2 shows ships which are in dangerous position 
because they are being acted on by radically different 
pressures. Reduction of speed increases the danger of 
collision. Also the pressure effects are more exaggera- 
ted and extra care is required in maneuvering (in depths 
less than 20 fathoms). 

It is understood that to maintain station during 


UNREP operation a certain amount of rudder is required 
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Wwhien depends on size, load, sea conditions, speed and 


ορ separation. But this change of rudder deflection 
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Figure 2. Ships alongside underway. 


decreases handling capabilities of the receiving ship. 

The problem of interaction effects has been studied 
theoretically by Silverstein [3] and experimentally by 
Newton [4] with both approaches showing agreement in the 
general conclusions. 

When two ships are underway with a certain separation 
and on parallel courses, the pressure fields mix which 
results in an unbalanced force and moment on each ship. 
This S moment must be cancelled by the rudder 
action in order to maintain station. However, there is 
a position where the rudder force tends to add to the 
force of attraction, see Figure 3. Thus in these posi- 
tions a bigger rudder deflection is required so the yaw 


angle produced in this way creates an outboard force 


TI 








Figure 3. Relative position where both interaction 
forces and moments tend to draw one ship 
toward the other. 


which balances both the rudder and interaction 
force. 

Newton's experiment, both with models and full scale 
trials in open sea, prove one important thing: It is the 
process of approach or departure from the abeam or "fuel- 
ing" position that includes the maneuvering risks. It 
should be noted here that an increase in velocity or de- 
crease of separation between ships results in the increase 


erfstcthe interaction effects. 
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ШЕ РЕКТУАТТОКЮ OF THE LINEAR MATHEMATICAL MODEL 


A. GENERAL CASE-SIX DEGREES OF FREEDOM 

It is well understood that the coordinate-axes system 
fixed on the ship's body is an Eulerian system. It is as- 
sumed that at t=0 the Eulerian system and the space co- 


ordinates system coincide, as shown in Figure 4. The geo- 
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Figure 4. Orientation of fixed axes and moving axes. 


graphic coordinate system then describes the motion of the 


Ship through the fluid medium in six degrees of freedom. 


ДЗ 





ШИС Newton's laws of motion in vector form are: 


— 
F (External Force) = ae (momentum) 
> d — — T TY ——  , ° 一 一 
M (External Moment)= dt (angular momentum) ( エエ エー] ) 


Then the six equations (six degrees of freedom) describing 


the ship's motion have been found [5] 

5 Ж? 5 . 
X=m [Ú+qw-rv=x¿ (5 tq )жу-(ра-г)%2о-(рг44)| 

1 2-2 Н : 
y-m[V*tru-pw-yc (r tp )+z.(qr-P)+x.(qP+%)] 

ξ DEMNM? š : 
2=m [w+pv—qu—z ¿ (Pp +q ) *x o (rp74) фу о (гр+р) ] 
=pI„+(I,-1,)qgr+ml[y (#—чи+ру)—2 „ (V=pwtrv)] 
ΙΤ νην). αατρν)! 


N-fI,t(I,-I.)pq*mixe(v—pwtru)-yg (à-rv*qw)] (III-2) 


Using (III-2) the following equations can be satisfied: 


zz 4 ñ ñ 
F = iX + jY + Ка 


(TEIS) 
= = Ë 1 
М = 1К + JM + kN 
and where: m — mass of the ship 
ОТТА Components Of force in the x,y,z directions 
K, M, N — components of applied moment about x,y,z 
axes 
а, V, м — components of velocity in the x,y,z 
directions 
Хе, % — distances of origin from center of gravity 


J my directions 





eee dyes components of angular velocity about the 
X,y,2 axes 


ΕΤ; — moments of inertia about the x,y,zZ axes 


Bauation (III—2) describes the reaction of the rigid 
meayeot the ship itself as a function of its geometric 
and physical characteristics. Note that they do not in- 


clude any external moments such as due to fins. 


B. HORIZONTAL PLANE MOTION 
It 1S well understood that the ship's motion in calm 


waters 1S described only by the following four equations: 


X = n [à*qu-rv-x, (x^«q^) «yc (pq-£) &z (peta) ], surge 
: 22 3 . 
Y = m[V+ru-pw-—-y (r +p )+z  (qr—-b)+x-(qp+Ë)], sway 
N = A (Y pwrru) y (ü rvtqw)], yaw 
K = БІ (1/1 jartm(y, (woqutpv)—z(v—pwtrv) |, roll (エエ エー4 ) 


because under the calm water assumption it is true that 
Roll = Pitch = Heave = 0. That is the horizontal motion 
Miles: Б = р = а = G = w = w = 0. Hence the equations 
(III-4) can be written: 


X = m [üu—rv—x r^-y t] 


G 
: . 2 
Y = m[V+ur+x_Y Yor ] 
K = m [一 2 (V+ エ V ) | 
A E UE エエ エー5 
N fI_+m[x saru) yg (ü rv)] ( ) 


Equations (III—5) can be further simplified assuming 


ӘШІР Ше center or Gravity (CG) is placed at the origin 
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of the x, y, Z coordinate system then γε. e = 0 
Hence, neglecting the roll equation we have: 

Xx = m[àü—rv], surge 

Ү = m[v*ur], sway 

N = ri , yaw (111—6) 


2, 

There are the reduced equations for steering and maneu— 
vering of a ship. It is noted that the left hand side of 
equations (III-6) represents the forces and moment along 
the coordinate axes, and the right hand side shows the cor- 
responding dynamic response terms on the horizontal plane 


Fer motion. 


C. LINEARIZATION THROUGH TAYLOR'S SERIES EXPANSION [5}, [6] 

The forces and ШІ on the left hand side of equa- 
tions (III-1) through (III-6) can be expressed as functions 
of properties of the body, properties of the fluid and 
motion. 

Since steering and maneuvering are of interest, forces 
and moments are also considered as functions of rudder 
{control surface) deflections and the change in r.p.m. 


(An) of the propeller shaft. Thus: 


Forces ) = tf (properties of motion, rudder deflection,etc.)= 
E. = f (Xo Yo’ Zo $r PU, VW, Pq, rÙ, V, W, D, QuE, 
Oe 
о епгаглоп MOCIOR рагатпесет 5 
parameters 


6,6,6,An, etc.) 
Ku R 
control surface parameters ( エエ エー/ ) 
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For a surface ship moving on the horizontal plane no 
Ir Ces or moments are due to orientation changes. Thus 
the forces and moments will be functions of three degrees 
of freedom of motion parameters, rudder deflection para- 


meters and changes in r.p.m. Hence: 


Y = f (u,v,r,u,v,r,6,6,6,An, etc.) (エエ エー8 ) 


In general we linearize a function f(x) by the use of 
Taylor's series expansion. 


n 
_ Ax af (x) Ax” a” f(x) 
f (x) = Ё (ху) zs PE NTC TE TF 1 n 


where Ax = x — х, 


For small values of Ax the second order terms can be 
neglected and thus considering only the following expres- 


Shon TOr f(X): 


Ly 
£(x)= £(x,) + Ax ατα 





GETI-9) 
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This is the linearized Taylor's series expansion. The 
same principle can be applied for small perturbations in 
Ener eauation (ІІІ-8), which is a function of many 
variables. 

Since the Taylor expansion is written for a particu- 
lar point we choose this point to be an equilibrium posi- 
tion. An equilibrium position is that of straight ahead 
motion, at constant speed with rudder amidships. The 
hydrodynamic forces and moment have been found [6] to be: 

x=XAurXv+tXr + Хайз Kev + үт + X 


5° + X An 


where Au =u == u; and ППП = ny and 

the subscript 1 is referred to the values of the variables 
at the initial equilibrium condition and where all tne 
partial derivatives are evaluated. For Y and N similar 
expressions hold. Equating the linearized expression for 


X, Y, and N with equation (III-6) results in the linear- 


ized equations of motion for steering and maneuvering: 


X AutX v+X r+X-UtFX-V+X-T+X.6+X An = mü 
u V га u END 6 n 

Y Au+Y v+Y r+Y-ú+Y-v+Y-r+Y.S+Y An = m(v+ru.) 
u V Те u V 7 6 n 1 

N Au+N V+N r+N-u+N-v+tN-r+N.64N An = I r AO 
u V s u V l^ Š n Z 


Note that the term mvr in the right hand side of the first 
of equations (III-6) has been dropped since the ship was 
assumed in straight ahead motion. The derivatives 

27. Lr at E 5 


Yue ns 
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and Nr Ne vanish for any symmetrical port and star- 
board shape of ship (symmetry about the xz-plane). This 
ἐπε effect of decoupling surge from sway and yaw. 


Thus equations (エエ エー10) become: 


(X-—m) u+X Au+X An = 0 
u u n 


(Y -—m) v+Y V+ (Y —_mmu,)r+Y.Xr+Y,6+Y An = 0 
V V iG Ë n 


1 Š 


(N-—I_)r+N v+N.v+N r+N,ó+N An = O (III—11) 
m" z V V r 6 n 


Using aS a basis the linear mathematical model of 
r ıons (III- 11) the criterion for dynamic stability in 


straight line motion [5] is evaluated as 
C = Y.N. = N (Y muy) > 0 


It is obvious that the treatment of the linear mathema— 
tical model, for the study of a dynamically stable ship, 
can be applied -to predict maneuvering-end control only for 
small deviations from the original straight line motion 


due to small rudder deflections and small changes in r.p.m. 


D. NON-DIMENSIONAL EQUATIONS 

Table I gives the dimensionalized, nondimensionalized 
quantities and their respective conversion factors. The 
prime notation corresponds to the nondimensionalized quan- 
tities. In equations (III-10)the force equations are 
2 


and the moment equation by (0/2) L7 ut. 


divided by EE 





The resulting nondimensional equations are: 


(X! —m’)ü’+X’A'+X'A'n= 0 
u u u n 


ЖЕП IU Y (У —m'u!l)r'+Y!r'+Y!16'+Y'A'n = 0 
V V r 1 r 6 n 


i 7 8 \} νι IE о с ес! ΤΑΔ na — s 
(N; I,)£ TNOV +N. V +N r +N 46 +N ὃ ПЕ 0 (IV—1) 


It is noted here that 


|. 
> 
el 
= 
© 


с 


for small perturbations. 
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Symbol 


TABLE I 


NONDIMENSIONAL NOMENCLATURE 


Non-Dimensional 
Form 


C = 


<l- 


N| = n|]. NS No] кој мк 


wir 


x 


DIE» 


pL u] 


pL u 


e N 
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Definition 


Derivative of longitudi- 
nal force component with 
respect to longitudinal 

acceleration component u 


Derivative of longitudi— 
nal force component with 
respect to longitudinal 
velocity component u 


Derivative of lateral 
force component with 
Respect to elansverse 
velocity component 


Derivative of lateral 
force component with res— 
pect to transverse 
acceleration component 


Derivative of lateral 
force component with res- 
pect to yaw angular 
velocity component 


Derivative of lateral 
force component with res- 
pect to yaw angular 
acceleration component 


Derivative of lateral 
force component with res- 
peet to rudder angle 
component 





e - 


He = 


No] Fe 


əк 


ο] 
で 
E 


tÑ] = 
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Derivative of yawing mo— 
ment component with res— 
pect to transverse 
velocity component 


Derivative of yawing mo— 
ment component with res 一 
pect to transverse 
velocity component 


Derivative of yawing com— 
ponent with respect to 
yaw angular velocity 
component 


Derivative of yawing mo 一 
ment with respect to yaw 
angular acceleration 
component 


Derivative of yaw moment 
component with respect 
to rudder angle component 


Derivative of longitudi- 
nal force component with 
respect to change in 
propeller rpm 


Derivative of lateral 
force change in with res- 
pect to propeller rpm 


Derivative of yawing com— 
ponent with respect to 
change in propeller rpm 


Yawing angular velocity 
component 


Yawing angular accelera- 
tion component 





мін 


Velocity СЕЗОЕТСЕІП Of 
body axes relative to 


fluid 


Transverse velocity compon— 
ent of origin of ship axes 


ге ве“ Бо Рака 
Transverse acceleration 


component of ship axes 
relative to fluid 


Hydrodynamic longitudinal 
force (positive direction 


forward) 


Hydrodynamıc lateral force 
(positive direction to 


starboard) 


Mass of the ship 


πο μα orr origin orf 
ship's axes along the 


x—axis 


Acceleration Of Origin of 
ship's axes along the 


x—axis 


Moment of inertia of the 
ship with respect to the 


z—axis 
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ΤΟ πι κ RELATIVE TO THE 
EARTH—LANGRANGIAN SYSTEM 
Newconvert equations (III—-6) from the Eulerian system 
to the Langrangian system, see Figure 4, the stationary 
Peordinate system, the total forces in the x and y direc- 
tion in the ship coordinates must be expressed in terms 


of the stationary coordinates as follows: 


> 
|| 


Xo COS U t Yo sin yY 


< 
|| 


x Е 
= ος Е sin ゆ (V—1) 


Also the center of gravity possesses velocity and ac- 
celeration with respect to the stationary coordinate sys- 


tem expressed as follows: 


x =u cos Jj) — v Sin Y 
O 
G 
у =u sin үу + v cos Y (V—2) 
O 
G 
End 
ae = ú cos y — v sin y —(u sin y + v cos V) uU 
G 
Y. = ú Sin Y + Ç cos Y +(u cos y — v sin ψ) ᾧ 
G 


S=) 


Не тола, (у“-2) сап be used to calculate the trajec- 
tory and velocity of the ship with respect to stationary 
coordinate system. Equations (V—3) can be used to calcu- 


late the acceleration along this trajectory. 


24 








V. ANALOG COMPUTER PROGRAMMING FOR 
ESCHTOHIP INDEPENDENTLY 
Pee PRELIMINARY STUDY 
In this section we will discuss the simulation of the 
amp dynamics in open calm sea. The Berne one Of motion 
for each ship are programmed independently and no interac— 
tion forces and moments are considered. 
1. Linearized Equations of Motion 
Equations (ІУ-1) аге repeated here for convenience 


as follows: 


(m'-X;)ü' = x u + X ón, Surge 

(m'—Y-)v' E RU +(Y mm use AY IS У бе , Sway 

(END = NV" + мн ну био "иу би, Yaw  (V—1) 
where 


r = ў 


т 1) 
ШЕШ oR 
Апо = оп 
πο = = 


It is noted here again that all the variables corres- 


pond to small perturbation variables. It is first chosen 
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е equilibrium condition as follows: 


Nominal values: Х = 15 knots = 25.32 Te sad = i 
у = 0 
O 
0 


2. Magnitude Scaling of the Linearized Equations 
Following the usual procedure each equation of 
equations (V-1) is scaled by letting each variable be re- 
presented by its maximum value times the scaled variable. 
That is: 
à = ü à 

and in non-dimensional form: 
' 
m 


where: Un is the maximum expected value 


Co 


is the scaled variable whose value ranges from 
gc Computer units (1.e,. + 100 volts for 
CI-5000 analog computer). | 

ЕЕЕ ЕН 5 substitution for all the variables in 


equations (V—1) and solving for the scaled variables ú, V, 


and 7 we obtain: 


E X, ur _ x 6n 2 
ШЕ ГГ! Dar ón 
ú” m 1 
Y! y! (Y !—m ' u фу ι i 
= Vv m Y O E i 
V Tn ΙΤ e u = 297 で ; 
V zu 
ME OP! Yo on 
ó R m ERG n Tan 1 
La R H aa 9n 
m V 
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хн | 
- N, V N We 


ф' = [ μα ο ραπ i 
τ МА MO 
NV N: ος. 
+ K | τι - ΕΠ 6R 
(1,-М ) T а 
N' On и 
+ (Y бп' , | (V—2) 
Еа h 


Letting each term in brackets be a coefficient (k), 


equations (V-2) can be written as: 


πι. -. .. “а ' 

uf = Куч + K」 っ 9n 

gc == TT — ーーー 

ў' = КУ Κα οὔ k, + Kig9R ko9n 

ҮГ — — m= umb pen 

Y Куу Κηγοῦ к1137 k4149R +K] 」=On (V-3) 


where the appropriate sign has been extracted so that each 
coefficient is now positive. Table IV gives the nondimen- 
Sional values necessary for equation (V-2) and (V-3). The 


following data:are also given in Table II. 


TABLE II 


с 
|| 


м2 ЕГ/Бес (nominal ship's velocity) 


O 
44 = 640.0 (£t/sec)‘; p = 2 lbs ΛΔ EX 527.97 rtt 
i32 _ | 2 _ 
ぅ PL ας = 9.428E+10; 501 u = 17.86E+07 
01 = 4.095E+13; 5017-1.4707Е407; 5012=2.78628+05 
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η. ο του οἱ Potentiometer Coefficients 
In obtaining a set of scaled equations, it is best 
to work with dimensional quantities. Actually, scaling the 
nondimensional equations is equivalent to scaling the di- 
mensional set. Tables III through V give all the necessary 
даға for the scaling procedure. Referring to dimensional 


quantities equations (V-2) can be rewritten as: 


E. X, Un X on, e 
u = I rm F beu on Surge 
ц m 
4 uv. (Y, mud 


ШЕШ туз !У * 19 


Ne Y z СТО ЕЕ 


F mv, ТА i TY ee J SR 


К оп 
Dr om 


ier. θη, Sway 
ν΄ πι 


τα A — Um No V. 


уп. Y 
Y πα. a ΠΛ ] V SE [1 Ex 3T. ] i aay)? 





NOR Rn T Ма Dm — 
+ Ши хуу !°^- + any on: Yaw (V—4) 
ze) m 72 ゆり m 
Let p = (m-X,)à. - mü. — X;à. 
= (1.166 x 10? slugs x 0.844 ft/sec^ 
21.0295 x ο. lieet sce x 0.844 TO SIS 
= 0.984 x 10° + 0.0866 x 10° = 1.07 x 10° lbs 
ру. - (m-Y.) v - mV. — ша 
208.0059 x 0.844 + 1.043 x 10° x 0.944 
= 0.984 x 10° + 0.880 x 10° = 1.864 x 10° ibs 
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da Y IN s к Uni NG 


ЖООШ ne fl/rad/sec- x 0.87 x T 


ШОО о КООШО; „ 10' — 2.777 x 10/ lbs-ft 


|| 


Then from equations (V—3) and (V-4) 





2 πι! у: [28.465 x ο τα. x 5.064 ft/sec 








K = [ ] 
11 D ШКО Ol ais 
БОО x 107 22 
Zee ZZ = 0.0623 x 10 - - 040623 
1.07 x 10 
X sn 5 rec/sec 
K БЕЛ Пп) .. [1:72 x l0 lbs/rev/sec x 30rpm x 60rpm ] 
12 Da O ps 
5 
Е 1.72 x 10 х 0.2 - 0.08037 
1.07 x 10 
NO 4 с 2 
K EI V ШЕ (9-16 xX ο ]5sZftct/sec x 5.064 ft/sec) 
15 D] っ 186 UN Le 
4 
но = 33.798 х 10° = -0.23798 
1.864 x 10 
ч (Y, mu) b, _ у те Un 
ュー | = (> 
ПО 12 
7 д 6-95 2 
B (1-042x10 lbs/rad/secx3.49x10 rad/sec-1.166x10 ft/sec x 
ОЕ pS 
O 3.49x10 “rad/sec, 
5 6 =2 5 5 
- [3.6365x10 —29.5x10 x3.49x10 1-13-6365х10 Ποιο | 
Ато” ООС S 
5 
БЕСІ = 5762 x 10 = —0.35762 
1.864x10 
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lbs lbs 








K . B " (72. 096x10^rad/sec^x0.87x10 ^ rad/sec^ 
17 D ШЕСІ. 
4 
ЕНСЕ) 09764х1072 - -0.009764 
1.864х 10 
B. ER m. М (4.823x10”1bs/rad x 0.35 rad, 
18 ЕТ? i о" 
© 
ИО 688 хо _ (090558 
1.864 x 10 
SA 4 lbs rev/sec 
K = n m _ ーー・ う 36x1 0 rev/sec x 30rpm 60 rpm 
13 > eo 
—1 936510: 0.5 --0 968х10“ 
= AA = ーーーーーーー = -0.005193 
1.864 x 10 1.864x10 
NS. Y 7 £ 
E = V m - [-1.307x10 fit^tiAceco5.064 ft/sec, 
ES 13 2.777 x 10! lbs-ft 
6 
Б 278 С 2.3834 ог -0.23834 into a gain of 10 
2.777x10 
Hio 。 lbs-ft 2 
K.. = [Um = [24.462x10 rad/sec x 3.49x10 rad/sec, 
ΤΡΙ: η, πο ΠΡ... 
7 
= Depo esx) = -5.6076 or —0.56076 into a gain of 10 
pu T0. 
Es | ; lbs-ft 5 
K...= [YM 2 (ШІ.552х10 ft/sec* x 0.844 ft/sec’) 
113 πε s 
2.777 x 10° 1bs-ft 
7 
1.30988х10 HN 
LO : 
lbs-ft 
ПОЕН I USSMIOP rad x 0.35 rad 
ENERO [ — στὰ 
їз 2.777 x 10° 1bs-ft 
-4.158x10/ | | 
NE = —1.497 or —0.1497 into a gain of 10 
PTT 1x10 
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N ón lbs-ft rev/sec 





m, . ,5.110x109rev/Sec x 30RPM 60RPM 
13 2.777 x 10° 1bs-ft 
6 6 
Ж БЕНО] x 10 o E ΠΡΟ ΧΤῸ = 0.0920 
2.777 x 10 ПП о 


The calculated potentiometer values are summarized in 
Table VI. 

Making use of the k's coefficients, equations 
(V-4) can be written as follows: 


Kk," + Κα 997 (V-4a) 


el 
|| 


kicV + Κι οὗ + kV + κι ОК On 


8 19 


y = K] 」 V+ к1127% k1134V* Кај ¿9R+k, y «δη 


Analog programming of these equations will give 
the left hand side quantities but sign inverted, because 
summer amplifiers and integrators invert the sign of 
their input quantities. Hence for analog programming 


equations (V—4a) are written as follows next: 


n= a u + к,» én] (V—4b) 
V = -İk s v + K16 ф + K] っ у + Кур QR + Kig én] 
y= -1К111У + Куто? + К 139 + Ky 449 + Κ11ςδη] 


Equations (V—4b) can be written as follows after 
taking care of the actual sign of each coefficient as 


already have been calculated previously: 


EU. з a —4 
ú Em u к,» ón | (V-4c) 
Ç = Kis V + куе Ú + Κ17 y 一 Кур SR + Куо n 

ђ = ә т SS < S 

ΙΙ ο E4459 * k4,13Y * K,,49R * K,,,5n 
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Equations (V—4c) are finally programmed in the 


Код С1—5000 сотпрцћег. 


CABLE ITI 
DERIVATIVE DIMENSIONAL VALUE 
7 1.042 x 10° 1bs/rad/sec 
7 2 
E —2.096 x 10 1bs/rad/sec 
Yy —8.76 x “ae lbs/ft/sec 
y 6 2 
ç —1.083 x 10” lbs/ft/sec 
N, -1.307 x 10° 1bs-ft/ft/sec 
N. ЕНЕ e 
Y. —1.936 x 10° lbs/rev/sec 
К: 10 2 
Ü —1.762 x 10 IDs-ft/rad/sec 
Yop 15623 X a lbs/rad 
Мер “-1.188 х 10° lbs-ft/rad 
X, —8.465 x ioe lbs/ft/sec 
5 2 

X. —]1.0295x 10” lbs/ft/sec 

X 5 

а 1.720 x 10" lbs/rev/sec 

Na ЛО: ο. lbs-ft/rev/sec 
B —4.462 x 107 lbs—ft/rad/sec 
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DERIVATIVE 


TABLE IV 


NON RD TUMENA FACTOR 


1 3 
z 2 


1 3 
ΖΡ; 
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NON—DIMEN. VALUE 


EE UE 


oe 


TOU T πο. 


η πη 
.0028 
ー.0027 
аз 
-0.0071 
-.0002 
BES 
-.0000052 
ー.00043 


. 0027 


=,00126 


= 0012 


25290007 


200/005 


—. 0000026 


—. 00227 


0079 





= 
|| 


= 
|| 


н 
|| 


ie 26 x 10° (slugs) 


TABLE V 


PARAMETERS AND MAXIMUM VALUES 


— mass of ship with entrained 


water 
m,/1.08 = 1.166(slugs) — mass of ship itself only 
1.43 x 10° Inertia (slugs — FT2) 
0.26 Yaw Angle (RAD) = 14.9 (degrees) 


0.0349 Yaw Ve1. (RAD/SEC)= 2 DEG/SEC 


0.00087 


5.064 


0.844 


550 


20° Rudder -Defelction (DEG) 


30 Propeller Speed 


160 


Yaw Accel. 


Ý 


om 


om 


527.85 Ship Length 


(RAD/SEC*)= 0.05 DEG/SEC* 


(PD SEC) 
(FT/SEC2) 


МЕСЕ) 


0.55. ιν 


(RPM) (RAD/SEC) = 


( REV/SEC ) 


NJ 


(ET) 


(ET) 
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TABLE VI 


ΠΕΙΤΕ VALUES 


Kii = 0.0404623 

Ki) = 0.08037 

Kis 20523798 

Ki 805095702 

K14 = 0.009764 

Kig — 0.090558 

Κο = 0.005193 

Kiii 2005594 (into a gain of 10) 
К 12 056076” into a gain of 10) 
X113 = 0.47169 

К 14 БИО (into a gain of 10) 
K] ] c = 0.0920 


For ship No 2 the following relationships hold: 


Em - K 


K っ っ = Κ12 etc. 
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TABLE VII 


MAKTNER CHARACTERISTICS 


Length, ft. 

Beam, ft. 

Marit, ft. 
Displacement, tons 
Ck Coefficient Ch 


Prismatic Coeff. E 


Midships Section Coeff. C 


528:3 
70.0 
29175 
207900 

56125 
. 6246 


.9807 


4. Analog Patching Configuration 


Bess Static Test 

MSRP NO” LT. 
patching used for the static 
chine was used for this job. 
tains the source deck of the 
Beant out. Table VIII gives 
used. 

(Ship No 2. 


patching used for the Static 


Figure 5 shows the analog 
eest ης 9300 digital ma 一 
Computer ae GO ロー 

program as well as the 


the values of parameters as 


Figure 6 shows the analog 


test. Computer program IB 


gives the source deck and Table IX gives the values of 


parameters as used. 


The static test was satisfactory since 


the printed values were the correct expected ones. As 
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| -199 


f д |- EJ Poly $. e) 
| ᾧ k | E 
D 1423 ne \ opa To6 i | 









Figure 5. Analog diagram for 
static test of 
ship No I. 


тег the coordinates transformation 
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TABEB VIII 


Potentiometer Adjusted Parameter Assigned I.C. 
Address Value of integrators 
FOr cta Iic test 

Ρ000 0217000 An A001 3. volts 
P001 0.0300 το ΑΡΗ ӘНЕ oO. volts 
P002 0.0803 к,» A005 9 volts 
P003 0.0600 IC—A003 A013 125 volts 
P004 0.0404 K); A0]7 15. volts 
P005 0.0052 Куо А051 3. volts 
Ρ007 0.0900 ТС—А00О5 A057 45.4 volts 
P010 00582379 K 

15 
P011 022383 K111 
P012 0.4717 Kr 
P013 0.0250 ри] 
Р014 0.1200 TC=A0T3 
P015 0.5607 K] ] 2 
Ρ016 0.1340 Г Р) 
Р017 0.1500 JC Ы 
P020 0.0920 К 15 
P021 0 "0097 K っ 
Ρ022 023576 Куб 
P053 0.0905 Кур 
P054 0.0500 
P055 0.1497 Kiya 
Ρ050 220337 I 
bol 0.3000 το A05 
P056 0.0092 D xu 
P057 0.4540 IC-A057 
P437 0.1000 AR 
P006 0.1660 Ио | 

п m 
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Figure 6. Analog diagram 
for static test 
of ship No 2. 


After the coordinates transformation 





TABLE IX 


Assigned I.C. to 
imeegrator for 





Potentiometer Adjusted Parameter 
Address Values 

Ρ024 0.1000 An 

Ρ025 0.0300 16-2035 
P026 0.0803 K っ っ 
Р027 0.0600 πο 02 
Ρ030 0.0404 K>1 
P045 020920 Κρις 
Р031 0.0052 K>o 
P023 0.0905 К.в 
Ρ034 0.1660 Га Ир 
P033 0.0900 LC—A033 
P046 0.4717 K っ ] っ 
Ρ032 0012379 K っ = 
Ρ035 022383 Κ211 
Ρ047 0.1497 K^14 
P036 0.0097 K^ 
P037 20.0250 εκ. 
Р040 01200 ТС—А041 
Р041 0.1340 [D Уу 
P042 0.1500 IC-A043 
Р043 0.5607 K っ ] っ 
PO44 0- 3576 Ko 
P407 ШЕМЕНІ BY) 
P425 0.3000 IC=A065 
P426 0.5000 

P052 055000 

P417 ο ου Be 
P427 0.5450 TC—A061 
P416 0.1000 AR 
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static test 


A035 3. volts 
A027 6. volts 
A033 DENN Volts 
A041 12. volts 
A043 jr volts 
A065 30. volts 
A061 54.5 volts 





Bar as the logic board is concerned Figure 6c shows the 
wiring necessary for the amplifier's operation where RT 
means real time and is connected Coa сала со О: ТЕ. 
This Figure was made for amplifier A001 and A003, but 
asis true for all the amplifiers used. Also, Figure 
6d shows the logic operation for the amplifiers A001 and 


20057 which is again true for all the amplifiers used. 


Апоов eat. £ Logie Goarol 
Br T4oo 


v- Io recorder 
Charmel 2, 
Όσο 
+100 yolks ο | 
AMS 1 | 


ENLCHARACTERISTIC LINEAR RESPONSE OP MARINER 

First of all, amplifier A011 was used as a timer for 
the κ. channel of the chart recorder. The wiring was as 
shown above. 

So by this configuration, a pulse was recorded every 
Meoee, 1.€., every one time unit. 

Figures 7 and 8 show the analog patching diagrams used 
to obtain the responses of ship 1 and 2, respectively. 
These response curves are shown in a later section. At 


this time interaction forces and moments are not included. 
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Figure 6c. 


R1 UO [- βοοι-ξ 
4. O &oo4-& 
DRWER | 
H4 PA A004-H 
Aoo- it 
oO.1KC 


s p 
τς 
wi) C4 
1O KC NAN) 
RT | 


Logic patching diagram for the amplifiers used. 


—O pooi-4.0 pr 
ar o— 
— (C Aoo ぅ ・ 4.0 pt 


A001 - 0.01 pF 
TS 
дооъ- 0-04 Р 


ARTS TS =D ZN REAL TIME 


={ IN МЕ SCALE 


Figure 6d. Schematic diagram of amplifier's logic opera- 
tion regarding the feedback capacitor. 
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Е F: 62:3 
| Steering (раб со 6. Ng 
P47 FSIK | 


Ee e 


Figure 7. Analog diagram for study of dynamic behavior 
Om snip No. 1. 
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| «63-69 | 
| 


| Ms 
rte > 


Figure 8. Analog diagram for study of dynamic behavior 
о пао Мо 2. 
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The dynamic behavior of the ships separately is of inter- 
Bete Also note that all the initial conditions were set 
ај to zero. 


fee Dynamic Test 


R C Amplifier A001—A026 
Cin р Co E 1 1 
с Сп" ве. то. 
$ s 
(ес! B SEE ΓΕ SES 
S "το 
A 1% 
ћ 1+ЕС S+1/RC 
1n s 
T= Bm l - l sec Settling time = T_ = 4 sec 
RC (1М9) (11F) 7 s 
e zem e 
Iri RC О 
s 
s κο 
RC tO S+o/ RC 
1 1 = 
= E wo. 94.8 
0.0404 0.0404 E 





(lo ΗΒ) 
Settling time = = 100 sec 


Amplifier A017—A043 


= i 
e == 2. 
O RC 
ME e р 
: т = κ αμ... 1 sec 
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Amplifier A013-A041 


R R R R 





R 
А, 1 T 1 1 E MC 
e€ = —[— e, + — e, + — e, + — e, + — p — ё |] 
O Ку 1 к, 2 Rio 3 Rio 4 Ry Rio O 
ë 
O 
O 
1ΜΩ 
= 0.1M2 
luF 
For the force free case it is 
е, = е. = 83 = e, Ξ 0 
Hence: 
е. = 108 és 
and -1 
ё = — /ов аж - 
O Кус O 
E мет і Ж 
=, μας g 10Bo "6 
1 1 
ΟΥ 1 | 
e T кшш с, let х=е. then х=е. 
у 1 
so x + ——(1080)x = 0 
R.C 
l 
Thus the characteristic equation becomes 
eo 
C 
1 
БЕ 1 RE 
πι 7 (Оу (0,5607) (0.0250) 14 Мне 
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Thus: 


at 
S 


By another way the same 


—é 
O 
ë 
О 
апа ес 
substituting 
ë 
О 
апа 
е 
О 
or 
Е 
(5) 
Непсе 
Е 
ШЕ) 
апа 
T 
T 
5 











Ен |- - 28 6 - 30 сес 
results can be obtained, i.e. 
= — — fe dt 
R) C О 
а a 
= — fe dt 
Ку O 
= — 1088, 
this back gives: 
Бе Q7 š 
= БС 108 f ĉo dt 
1 
OS : 
RC ) “о at 
JE 
ο о as) 
RJ C s 
SR,C : 
ЗРНА 10αβ _ SR」C+10oB _ 1008 + 1 
RC. SR,C SR,C 
10αβ 
10αβ O) (0.5607) (0:0250) 0.14 ` 
= 4x 7.15 = 30 sec 
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Amplifier A005-A033 
КООШО ОТ ene force free case only the closed loop 


is examined. 





x = =] Ку wdt 

x = KW 

w= —[f + k,10x] 
em for f=0: w = —[k,10x] 


Or Ба = κιν = —k, [k,10x] 


x + k,k,10x = 0 


The characteristic equation is: S+10k,k, = 0 


I2 ШЕН то ποιο IO sec 


ЕЗ 
II 


са 5 = 10 sec š 


Flgure 9 shows next the response when approximately 
ПИО спссек angle is applied. Note that + l computer 
GOmGS5omdS Eo + 20°. So a pot setting of P 437 for 
lini] and P416 for ship 2 equal to 1.0 corresponds to 20°. 
The sign of AR will be determined of course from the sign of 
the reference voltage. In fact for every curve taken by 


the chart the following relationship holds: 


+ 100 И а номе x sivolts/divisions] 
x lin Gna ts | x[in units] 
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Figure 9b. Dynamic test of Mariner. 
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Figure 10b. 
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where: = maximum expected value in actual units used 


in magnitude scaling the equation 


A = divisions measured 

SEU sensitivity in volts/division 

AS actual value of parameter x in actual units 
100 volts = 1 computer unit 


In Figure 10 the time scaling is changed and the 


linear response of Mariner is obtained in terms of u,V,V," 
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Figure 12. Steering Control used for the dynamic test 

of Mariner. 
when AR = 15° is applied. In Figure 11 the response for 
ПИО (5 75поуп. The time constants estimated pre- 
viously are the same as the obtained ones. | 

2. Steering Control 
a. Analysis of Steering Control 

Up till now a simplified model for the repre- 
sentation of the steering control was used, which produces 
τρ commands 4 or — according to what reference 
voltage is chosen. The steering control configuration is 
shown in Figure 12. 

' One question which arises immediately after 

looking at Figure 12: is this a physically realizable 


steering control? The answer is no for several reasons. 


SO 





Primarily the actual steering has a certain time lag in 
its response. Also for a given rudder change angle by the 


Eumsman the actual rudder accelerates at the beginning 





ἣν". η 
Gain of 4 


= 104 


and then decelerates at the end of its final position. So 
to be closer to the physical situation of a steering con- 
trol a definite time lag for the rudder action is needed 
aS well as certain limits on the rate of change of the 
mueder angle, but without changing the initial gain set. 
The time lag can be produced by the above shown circuitry. 
εὐ back gain is: k = 0.175 x 0.286 = 0.05 


The time constant and settling time are next calculated 









о/о EN 
O RC in RC О 
Е S 
ео e ο LORC 
6. s+k/10RC 
in 
-k = 0.5 
t = 1/(k/10RC] = 2 sec ς 
P 4 
= = 6 Initial slope 
m 4x2 8 sec = г 
Steering 
στ τὴ (ον οί 
Ri ER, 


“8 seu t 


Now regarding the above shown diagram for the rudder the 
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following can be deduced. The desired time lag is obtained 
as well as the big initial slope of AR. The desired res- 
ЕСЕ іс the dotted line curve AR, seen in the previous 
sketch. To achieve that a limiter 1S used which limits 


the output e, as it can be seen from Figure 13a. 





DEN >lo R 





Escure 13a. Improved steering control. 


Bhesrunectron of the limiter is similar to the 
Saturation effect as next shown. Initially the error sig- 
nal [$9e/6e, ] mI. ene which Corresponds to +100 
or -100 reference voltage, since the feedback is zero. 
Recht after that due to the feedback the error signal is 
reduced almost linearly to zero. The time required for 
this is the desired time lag. Figure 13b illustrates the 
various steps made for the improvement of steering control 
in terms of the output rudder command. Figure 14 shows 
the complete analog configuration of the steering control 


system as used for both ship models. ος 
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2. Command with time lag 
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3. Command with both time lag and rate of change 
limitation 


Figure 13b. Improvement steps of rudder commands. 
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b. Time Lag and Settling Time of Improved 
Steering Control 


From Figure 13a it can be seen that: 


О – а е; dt 

= E [еде 

EI Ете. a 10 out E 10 біп where 
lm 
INE EU S c Sg 
im 

Be -Ῥ E 

“out T5 E 10 out x 10 SE dt 


For the force free case and within the specified limits 





aM 15: 
e ニー Ја е a 
out ITO σας 
ë E e 
out JO OUt 
or ë a е — 
Out OS tout 
Miereharacteristieoc equation is: 
ap _ 
S + 10 0 
and て =1| [op/10] 
Now for: 6R_= 20° 
m 
de = 7? 
k = 0.5 
B = 1 for real time simulation and 
SR 
from E = m Е 
m 
Кбе, 
а = Ta ОИ 50 
m 
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2 sec 


mers obtained: T 


ll: 


and а 8 sec as before. 


Еге 15а апа 155 show the responses for input commands 
MEN und An = + 9 RPM but with different time scaling 
respectively. Figure 16a and 16b show the responses for 
mot commands AR = + 10° and An = + 9 RPM but again with 
w erent time scaling of the chart recorder. It is ob- 
served that the expected time constants and settling times 


can be measured. 


TRANSFORMATION OF COORDINATES 
MN Tntroductory Discussion 

The responses observed in the previous sections 
were taken with respect to the Eulerian system, i.e. moving 
coordinate axes on the ship. This makes difficult the 
measurement of the separation distance between the two 
ships during ene UNREP operation. It is desired therefore 
Lo use the space coordinate system, which for ease of com- 
puter print out is made to move ahead at 15 knots. It is 
also noted that the ship initially moves straight ahead at 
75 knots. Figure 17 illustrates the above discussion. 

ποπ спе inspection of Figure 17, equations (V-2) 


can be recalled, i.e. 


x 


5 DI τ υπ 


y u sin y + u cos y 
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Figure 15b-2. 


Dynamic test of Mariner. 


An = Шш 9 RPM. 
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СШ ЛАС point computer unit XDS-9300 is used in hy- 
brid operation, needed for the coordinates transformation, 
meen the analog unit CI-5000. Computer program IVA and IVB 
show the source decks used and Figure 18a shows the control 
statements used for both digital and analog computers. 
Figure 18b shows the communications trunk lines used between 
digital and analog computers. Figures 19a and 19b show the 
analog patching diagrams used for the coordinates transfor- 
mation for both ships respectively. Table X contains the 
potentiometer values used. It should be mentioned here that 


Meche control of the hybrid operation for ship No 1 the 


logic switches "0" and "1" were used. For ship No 2 the 
τ Switches "2" and "3" were used. The hybrid operation 
A RA 


ee == 2m2 


— P. 





A 


Figure 17. Space coordinates system 


was performed with both the switches being in the middle 


Position. 


70 





É 






TEST(Z) у» 444 


m ἜΝ fr En AA 


ту 
Í 
Figure 18a. 


ГЕКЕ ЕЛА сяго controllof analog and digital 
Wie 





ADC DIGITAL OPERATION PAC 


(9) 7590 ++ ship Wed 
A) 154 0, 


%, u St No 9 (Се, 
(2) 1502 Y, 










< 
t 


А ανα (2) 


у 


ха = Wes ΤΌ γε οφ 





πη 
Yo, Др (5) 


O Tek ie 7 
(5) ote oe 7o, oan 
(6) T506 m. 


%。 こい 。Cos 9-9, SIN. x. mt SNE) 
2 
Figure 18b. Communication trunk lines between digital and 


analog computer used for the coordinates 
transTrormatıon. 
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Figure 19a. Analog diagram used for hybrid simulation 
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Figure 19b. Analog diagram used for hybrid simulation 
for ship O 2. 
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2. Computed Linear Response of Mariner after 
@Gordinates Transformation 





Figure 20 shows the linear response of the Mariner 
 ртохтпасеіу + 10° changes referred to midships in the 
rudder angle. 

Initially the Mariner is moving straight ahead at 
ENEnots with the rudder at midships (6=0), and then 
EN UO" 1s applied for 25 seconds. Three phases are 
Observed for a turn. 

a. Initial Approach Phase in which 

V-0, yr=f=0, v=0, r= )=0, Y=0 
be First Phase in which 
r>0, і-с-)<0, v=0, r= p=0 
This phase Starts at the instant the rudder 
is turned, and ends before the rudder reaches its full 
deflection angle. The rudder force (8) and rudder mo 一 
ment (N Š) are dominant and produce accelerations which 
are opposed only by the inertia of the Mariner. 
' 7; ποισα here that the rudder force 
Y,9 15 towards starboard, since the rudder is at the 
Сп OL the Mariner, which corresponds to a port (nega- 
tive) turn although the transverse velocity v is positive 
(towards starboard). 
c. Second Phase in which 
v>0, v>0, к=фр<0, r=)<0 
It must be mentioned here that a force Y. V 


EN cated towards the center of the turn (towards port) 
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Bigure 20-2. 


Characteristic linear response of Mariner 
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which is balanced eventually with the outwards centrifu- 
Orce of the Mariner. The acceleration Y is zero 
when this balance is accomplished. 

The third phase of the turn is not developed 
because the rudder is given a negative —10? angle and 
the Mariner starts to turn Starboard. It should also be 
noted here that at the instant the rudder is at ~10° the 
center of the gravity of Mariner is at maximum port (πθσᾶ- 
=) position. The process of turning +10° and —10° 
angle to the rudder was performed several times. 

Figure 21 shows the linear response of the 
fee ue to + 5° angle of the rudder. The difference 
ce response obtained with + 10° rudder is basically 
the smaller amplitudes. 

Figure 22 shows the linear response of the 
Mariner due to a positive change of R.P.M. of the pro- 
peller speed of approximately 9 R.P.M. with respect to 
puesegurilibrium R.P.M. (corresponding to 15 cU 

The reference axes PS nen are moving at a 
Louoerty 15 knots initially Straight ahead as is the 
Mariner and consequently the (x,y) axes. The response 
is recorded by curves of calculated changes in the para- 
meters ф,ф,ф,у У 'Хо' Хо“ πο не ја те origin Of the 
(x,y) axes in the Mariner with respect to the(x Yo) axes 
250 ft on the x, axes and 48 ft on the Yo axes. Note 
B rexrhat the computer unit, i.e. 100 volts, equals 160 ft 


for yg axes аза Lor x, axes respectively. 
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Figure 21-2. 


Characteristic linear response of Mariner at 


Мен 5° 
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Figure 22-2. 


* 9 RPM 


Characteristic linear response of Mariner at 


An 


8'3 





The movement in the direction Χο is due to the X An term 
 ИСапасіоп (V—1). The origin of the (x,y) axes in the 
Mariner moves to the origin of the Cie) axes in 110 sec 
Bpproximately. The ship also yaws —1.5? approximately. 

τ ΤΕ посеа here that only the force 
moment due to the RPM propeller change are simulated. 
The force and moment due to the m or equilibrium 
speed of 15 knots is neglected. 

Figure 23 shows the characteristic linear res- 
ponse of Mariner obtained in terms of the calculated para- 
meters An,xo,X ,ў ,у ' (а У ап οσα dr OPM echange of 
Eu Tier speed. Again the Mariner is assumed moving at 
U S at equilibrium position. The initial position of 
the origin of the (x,y) axes was the same. 

Figure 24 shows the linear response of the 
πο to a rudder change of approximately 40.562? (for 
P437=+0.0281) and a simultaneous +9RPM propeller speed 
change in terms of the calculated parameter perturbations 


КИ х, У 


E X. Eyo сатав initially the Mariner 


ο 
is moving straight ahead at 15 knots with rudder at mid- 
ships. Then both the rudder and propeller speed change at 
+0.562° and +9 RPM respectively. This results in yawing 
moments which cancel each other (y=0). The only obtained 
eee bation 1S a change of location of the origin of the 


(x,y) axes on the Mariner but only along the x, axes. 


Thus it is concluded that a rudder angle of about +0.562° 
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Figure 23-2. 
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Characteristic linear response of Mariner at 
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is enough to create a yawing moment necessary to cancel 
ШОО lng moment created by a change of +9 RPM of pro 一 
ЖИЕ ЕЕЕ среейі. Figure 25 shows the response of the 
Mariner nea in terms of the calculated perturbations 


ο ο л, ДУ, АБ апа 


‚х,у О 


of the parameters An, x n 


O O 


created by a simultaneous change in rudder angle and 

BEEN Uer Speed at + 5° and + 9 RPM respectively. These 
Curves can be analyzed analogously with the preceding 
curves. It is mentioned here that the same response was 
K q for both ship No 1 and ship No 2 in all the pre- 


VlOus cases. 
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Figure 25—2. 
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INTERACTION FORCES AND MOMENTS ARE INCLUDED 


ENT RODUCTORY DISCUSSION 

During the UNREP operation the tracking (receiving) 
Ship is usually keeping course to avoid collision and 
maintaining station relative to the leading (replenish- 
ing) ship. The leading ship must keep a steady course 
and keep oscillations about this course to a minimum. 
Meretnis reason the remainder of this study is subdivided 
into the following sections, namely: 

1. Phase I: Interaction forces and moments are 
lieu oq ила аке applied only to the 
leading ship as it 15 overtaken by 
the tracking ship. 

2, "phase II: This phase considers the complete 
hybrid simulation of the UNREP 
Operation inciuding all the inter- 
action forces and moments. 

ШИ Осо улагыла phase 1 the response of the lead- 
ing ship is of interest. Further work in extending phase 
i will be the response Ce rien cj ig Ship'as 3t over- 
takes the leading ship so interaction forces and moments 
are applied Bene әтір. ΠΠῚΞ will be better 
puserstood after reading of next section, since it is 
related with how the interaction data curves are given. 
Briefly it can be mentioned that these data correspond 
το that case in which the interaction force and moment 


ESnoshrp A as ship B overtakes it. So as a first case 


Ship A is chosen to be the leading ship. 
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B. INTERACTION DATA CURVES-INTERPOLATION PROGRAMMING 
Figures 26 and 27 show the nondimensional N moment and 
Y force respectively versus the longitudinal separation 
(a) between midships and the lateral separation (f). The 
@emeratioOn of these curves by the analog computer exceeds 
its capability. Therefore these curves are made piecewise 
linear and stored on a digital computer namely the XDS— 
9300 machine. It should be mentioned here that in this 
K array the points which correspond to 600 ft longi- 
tudinal distance (a) are included for both N-moment and 
Y-force. However these points are not shown on Figures 
26 and 27 respectively. Computer program III shows the 
program necessary for a two-dimensional interpolation in 
two separate arrays for N-moment and Y—force respectively. 
It is mentioned that the longitudinal separation distance 
(a) can give negative or positive values but the lateral 


distance between ships (8) is restricted to positive 


- 
- 


values only. Figure 28 shows the geometry as far as 
Pongitudinal and lateral separation distances are con— 
S It should also be mentioned that in phase I 
interaction moments and forces are applied only to 
ship A (see Figures 26, 27, 28) which corresponds to 
the leading ship because the data curves for the 
Buberaction effects are such. To apply interaction 
moments and forces to ship B, which corresponds to the 
tracking ship, the longitudinal separation distance 
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Figure 26. Dimensionless N Moment vs. Longitudinal Separa- 
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Eure 29. Geometry of longitudinal and lateral separa- 
tion distances. 
appropriate value for M and Y respectively. M and Y 


are vectors with opposite direction in the case of ship 
B as compared with those applied to ship A. Ships A 
and B are defined in Figures 26, 27, 28. 
C. MEASUREMENT OF LONGITUDINAL AND LATERAL 
SEPARATION DISTANCE 
1. Longitudinal Distance 
To measure the longitudinal distance, a, between 


midships the following equation must be used: 


ο ο ος 
1 2 
oj ες: (VI—1) 
2 να. 
СОБЕ ЕЕ па in this equation 
а = а а 
m 
SZ Xo Xo 
2 2 2m 
zu, x Xo 
i 1 lm 
where x d уа are the maximum expected values of 
O O m 
2m lm _ T С 
these variables respectively, and XQ Χο εᾱ are the 
| | 2 
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Scaled variables whose values range from —1 to +l. Hence 


ο ντους (VI-1) become: 








—а, αι = x, X, E 25 (VI—2) 
2 2m TE Tim 
or x x 
οτι τς A y (VI-3) 
O 
m 2 m 1 
Choosing: ал” ο τε 
x — 550 It 
O 
2m 
x =) 550 4h t 
O 
lm 


equation (VI—2) becomes: 


— = x pr (VI—4) 
2 1) 


Figure 29 shows the analog patching diagram necessary for 


equation (М1—4) 


FROM A057 Xos 
i Digitol input 


“πσο5- ШЕ) 





— — FROM A064 


Figure 29. Longitudinal distance analog diagram. 


2 Lateral Distance 
Memmeasute the lateral side—to—side distance be- 


tween ships the following equation must be uscd: 
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EZ SAT 
Or —B = y ο (VI—5) 
ој O, 
Substituting in this equation: 
Бар В 
Y, = у у 
ke Nm 
y. = y. y 
=> の の 
where: y N ‚ ß_ are the maximum expected values of 
Отт ?2m E 
these variables respectively, and Yo ; Yo EC 
1 2 
scaled variables whose values range from —l to +1. Hence 
equation (VI-4) can be written as: 
1802-09 Y Y (VI—6) 
От 91 E» 
Or D y 
-B = [2 y, - L3 y, (vI-7) 
m 1 m 2 
Choosing: В = 100 ft 
С БО ЕЕ 
| 2та 
nm = 160 ft 
lm 
equation (VI-7) can be written as: 
-B = A — ΙΙ, (УІ--8) 
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Figure 30. Lateral distance analog diagram. 


93 





Eure 30 shows the analog patching diagram needed for 
equation (VI-8). 

It should be noted at this point that because 
me coefficients of Yo, and Yo, in equation (VI-8) are 


EE слап l, the pots are set as follows: 


P040 


νυ ποπ“ Gain of 10 
P052 = 0.1600 into a gain of 10 

During the simulation 8 can take negative values 
but the data curves are given for positive values of B 
only. Thus in some way an absolute value should be pro- 
Medea either by digital subroutine ABS or by analog-patch- 
img Configuration. Тһе latter is chosen since the use of 
digital machine adds time delay in the simulation which 


Maye give rise to an error roughly of 5%. 


—À 


FROM 4056 ES. 


— 0; いて し inp Л 


FROM fose Py 75.20 = (20) 


- 





Ше 314: Absolute value patching diagram of |]. 


Figure 3la shows the analog patching diagram neces- 
БЕКУ to obtain the absolute value of B 

ШЕ should be noted that in Figure 3lain the second 
ШИРЕ ОБ А024 summer the gain is 10 but since: P405 = 0.2000 
1η becomes finally equal to 2. Also the limiter L 10 
limits the output of A066 to negative values. More speci- 


fically —|8| corresponds to +1 computer unit limit and 
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+ ЕЙ ил соропаз to —1 computer unit limit respectively. 


However according to the following illustrated fact: 


N 





~p Me K+ NERO 


there is a better way of getting a gain of 2 since two 
resistors of value 1.0 M in parallel gives a value of 0.5 
as value of input resistor which corresponds to a gain of 


2 as shown next in Figure 31Б. 


FROM #056 =)” s 


~ Pipital input 


(3 
` EN 4 N xag Гы ~ ⁄ 
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Figure 31b. Absolute value patching diagram of В. 


Finally the way of Figure 3lb is chosen being more compact. 
EXUAPPLELICATION OF INTERACTION MOMENT N AND FORCE Y 

BOSSHIP A (LEADING SHIP) 

Ix Itroductory Discussion 

Ship No 2 is chosen arbitrarily to be the ship A, 

that is the leading ship, on which moments and forces due 
BOminteraction effects would be applied. Figure 32 illus- 
trates the way N moment and Y force are interpreted by 
hybrid operation of the machines. 

The digital program which takes as data the values of 


8, a, in the computer mode of hybrid operation, finds the 
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K I DITO UNTI  xPE- 9306 
ADK INTERPOLATION PROGRAM DAC 
FROM 1503 ‘Yo T425- (+) 
τ.) 
КБШ 2-77 


To T427-(8) 





FROM "T 507 To T422-(3) 
(7) 
Нисе зЗ2а. Interpretation of N-moment and Y~force for 


hybrid computer operation. 


corresponded values of N and Y and then gives back to the 
analog the values of N and Y, which are scaled quantities, 
Lakes care of: 
a. Rescaling of 8, a since B = 1.0 апа а = 1.0 
corresponds TCO 110 ft and 550 ft respectively. 
b. Making linear interpolation in the two- 
dimensional array. So it finds the corres- 


ponded nondimensional quantities of N' and Y'. 


G Tuncer rom the data curves, Figures 26, 27, 
N' = 44.5 x 10> 
max 
y' = 90.0 x 107? 
max 
hence, finding the dimensional values of N 
η Y namely 
max 
PE 15554 - ал 5 ~ 10 9.428% 107 
max max 2 O 
y 


ο LO LDS 一 も 
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AA 


ーー ニテ 『 ーー 
and Sax > тах s a ч 
ТОКЕ 737787 5 10/-1.608 x 10?1bs 
= 72.9 tons 
а. Finding from N' and Y' the dimensional values 


of N' and Y' respectively as follows: 


UN oL uf ns 07 
MY x αρα” ο ο οκ Don 


ERU Scaling both N and Y quantities taking in 
consideration the fact that the subroutine 
DAC multiplies each number given to the analog 
ЕПЛ БУ O Sí actor ot 100 1.e. it converts 
τη feonmputer units, so: 


NSCALE 


|| 


it 


YSCALE 


ах 


Computer program VA contains the source deck as 
well as the potentiometer values and amplifier addresses 
used for phase I. 

Computer program VB contains the source deck used 
An phase I, but ship A (leading ship) was in that case 
the second model ЕМ лү лет The difference was only 
in the trunk lines which were used in the previouse case. 
ШОШО 32 the trunk lines T426, T427 belong to the ship 
model No 1. 


2. Dynamical Representation of Interaction Moment 
and Force-Analog Programming 


Recalling Newton's laws of motion, equation (III-1) 


ШОШ Се the study is done for the ship to be in 3 degrees 
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Meeerecdom, equations (III-6), (У—1), (V—4), the applica- 
tion of the external interaction N-moment and Y-force are 
going to be an additional term in the yaw and sway equa 一 
tion respectively. Considering the dimensional form of 
the linearized equations of motion, equations (V—4) will 


become for ship No 2 




















= um = Хх бп ーー 
u = E ur [ Jp Pour gë 
D 2 
— NV ο л җи у ος NW u 
ШИ Y LL e 6g 
p 22 22 
Ἱ OR πο τομ 5 
+ E ÔR + ο ön + , Sway 
22 22 22 
[ЕКИ ο 
additional term 
л — Ка の NM EX 
0 = 10-0) у + [50 фо 10-9] v 
23 23 23 
DM OE Ол 
| R Po) т Ln , Yaw 
D D D 
28 23 n 
addicional term 
( V エ ー9 ) 
where: D era ] 07 x Ten 
в DI u m 
= n. 6 
D^» — (m ur S Ξε; "τ IDS 
Е "NM 7 ー 
D っ っ = (I, ο... 2.1777 x 10 1bs-ft 
Be interaction force 
N = interaction moment 


It is clear from equation (VI-9) that the surge 


equation remains unchanged to interaction effects and 


ο 





only sway and yaw equations have cach of them one addi- 
па! term respectively. Thus for the analog program- 
[emg OL equation (VI-9) it is only necessary to add the 
ENcusitry Shown in Figure 33a to the patching diagram 


лаике 19b. Note in this Figure: 


23 2777x109 


B here is a factor of division by 10 since the value 
(N ax! P23? > l and so it is necessary to feed amplifier 
Ern a gain of 10. Note that since in amplifier A036 
we need one more, not available, input gain of 10, a 


Separate internal resistor of 0.1 MQ is used through the 


Summing Junction (SJ) of A036 giving the required gain of 








ШЕ Also: 
τν; = s= = 1.608x10 _ (PEE QJ 
22 10172170 
ΠῚ 
Te verorder ΑΝ З 
42 | | 
м 
б МГ й 
T ed 5 
» M Er “09 = 
Е > μυ. T +25 (3) 
E 4036. EM o | 
(gain od 10) 
= a Ka | 
lo recorder 
ἘΠ A034- a - A207 A. Digital output 
(gxin of 7) ы. NS る (+) 
б; 
Figure 33a. Analog Diagram programming of interaction 


N-moment and Y-force for ship A (No 2). 


104 





For the ship model No 1 the analog patching cir- 
curtry of Figure 33b was added to this of Figure 19a 
whenever it was desired to be used as leading ship and 


ship model No 1 to be used as tracking ship respectively. 


pL M | Al 


го 
ο N E 
40714 
ET. ES 1016 , 


T4286 - (7 
(gon of ος BE ; C7 


= Dipital output 





Το recor A 16) 


Into 4006 а aaa zy ela ne, АТЫ 7 Li BU outipux 


Das т 4-2 p (8) 





Figure 33b. Analog diagram programming of interaction 
Dr onenesandey-torce Tor shrp model No 1. 


It must be mentioned here that for Phase II where 
interaction effects have to be applied on both ships, both 
the configurations of Figures 33a and 33b must be used. 

3. Obtained Responses for Phase I : 


M CaCl onary runs (1.e. both ships have same 
propeller speed) 


Figure 34a,b shows the obtained linear res- 
Fonse Of the leading ship (ship A) in terms of the calcu- 
lated parameters Y, N, Yo” a, B, v, V, ф 
Originally the leading ship was placed at the 


origin of the space coordinate system ος) and the 


tracking ship was placed at a 15241 Et and E 170 fr. 
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Linear response of leading ship for Phase I. 


Broure 34b. 


+524 ft, p= JO.ft) 


(АВ-0, ӧп=0, а= 
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Then both the ships were left to run with the same propel- 
ler speed which was 15 knots. The same speed of 15 knots 
was assumed for the origin of the space coordinate system 
Ey. 

The leading ship was forced to yaw to positive 
angle due to the additional term of the interaction 
moment N. The interaction force Y produced a lateral ac- 
celeration and lateral velocity which forced: the leading 
ship to move laterally decreasing continuously the lateral 
separation distance. In the sway equation (VI—-9) not only 


the Y—force term was acting but the derivatives Doux - 


p y 


were producing additives effects since due to the yaw 
angle a change of yaw angle and rate of yaw angle change 
were obtained. It is also noted that the interaction force 
Y and moment N are slightly increased during the run since 
the lateral separation was decreasing. Similar results 
were obtained for the second run in which the leading ship 
was again placed originally at the origin of (x, Ly) axes 
and the tracking at I ウン (hua Mian a +70 ft. Again 
Buc Origin of пе (х,у) was moving at 15 knots as well as 
both the leading- and tracking ship. The response was 
calculated in terms of the perturbed parameters, %,у,М,0, 
Bzv,v as shown in Figure 35. The lateral distance during 
this second run was decreasing faster since the values 

for interaction force-Y and moment-N were greater. Ob 一 


viously from the data curves both force Y and moment N 


EXconot symmetrical about the origin. 
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Figure 36 shows the response of the leading 
ship when the longitudinal distance and the lateral dis- 
пес were set equal to zero and 70 feet respectively at 
the beginning of the run. Now because the longitudinal 
distance is zero, the Y-force and N-moment have not only 
greater values than previously but the leading ship yaws 
to negative yaw angles. This has as effect a negative 
mal velocity and an acceleration initially which 
eventually reduces the interaction force-Y and moment-N 
drastically, but the lateral distance is increasing now 
Since the velocity along the y Taxis has opposite sign. 

For all the above three responses the rudder 
EN oh the Ships was set at zero, i.e. = 0. It should 
σος at this point that the hybrid simulation program 
(computer program V) was made to work between lateral 
Separation distances 100 ft and 50 ft since data are given 
for this range. The increase of the lateral distance is 
ame: 

a. To the negative yaw angle. From equation (V-2) 

it seems that eventually the term (u sin y) 

becomes predominant over the term (u cos 1) 

hence even if a positive value is obtained for 

the v velocity a negative value of y_ is 


actually applied to the leading ship: 


b. The negative yaw angle is obtained because for 
0-0, В=70 ft a negative moment N is applied. 


It has been mentioned previously that the two 
Models are identical. So for Phase I No.2 model has chosen 


to be ship A, i.e. leading ship. To complete our work for 
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Phase I and check the wiring for model No l the case was 
reversed and ship A was the No 1 model and ship B the No 2 
model. Remember ship B ls always the tracking ship. For 
this last case a change in wiring on the analog board was 
made for the circuit representing the absolute value of 
the lateral separation distance (ß). The change was due 
to the lack of operation of the limiter L 10 as was shown 
Figure 3la and 31b. 

The following configuration is used from now 
on for the absolute value of the lateral distance (ß) as 


is next shown in Figure 37. 


FROM 4056 


»- [o digital input 
ЗОИ) 






where De are solid. 


stale oL iodat g 


Figure 37. Absolute value patching diagram of the lateral 
distance. 


Just briefly is mentioned that it is desired 
бошсеер the voltage of point A between the limits 0 and —. 
This is because diodes D, D, are connected in such a way. 
Now when point C has a positive voltage D; is OFF and 


point B has a negative voltage hence point A has a positive 


voltage since D, Ondu ting. When point C is negative, 


Ip? 





voltage D, is ON and B has about 0 or +0.5 volts voltage 


a 


much is not enough to break down D, hence point A is at 


2 
zero voltage. 

Figures 38 and 39 show the obtained response 
for the leading ship, when the longitudinal distance is 
+400 feet and —400 feet and the lateral distance +70 feet 
from the tracking ship. The response was obtained in 
terms of the perturbated parameters U, y, N, "s απο, 

у, У . The ships were moving at the same speed of 15 knots 
ШЕП the leading ship placed at the origin of the ο ας 
ЕЕЕ Іп Figure 38 the tracking ship position was at an 
£400 ft and n +70. The interaction force -Y was almost 
zero during the time of this run, while the interaction 
moment —N had a positive value producing a positive yaw 
angle. Because ος was positive the track of the leading 
ship reduced the lateral distance (ß). 

ner discussion can be made for the response 
of Figure 39 where the tracking ship was placed at Ба 
—400 feet and Ne 70 feet. The lateral distance was again 
reduced continuously. 

Several "stationary" Buns а = ЗЕ сре same 
speed of 15 knots for both the ships, were made for the 
ngi tudinal distance (a) being 160 ft, 300 ft, 200 ft and 
the lateral distance (8) 70 ft, respectively. Figures 40, 


41, 42, 43, 44, 45 show these responses in terms of the 


perturbated parameters of the leading ship, WYN, Yor ar By 
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Figure 41-2. Linear response of the leading ship for 


Phase I. ο ου ου Зас асое ΒΞ +70ft) 
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v,v. The analysis of these response curves is analogous 
to the previous ones. The lateral distance in some cases 
is reduced and others increased by the interaction effects 
depending on the relative longitudinal position of the 
ships since the lateral distances were kept constant at 

ου τι at the beginning of each run. 

τοσα Еле Observation of these "stationary" runs 
it seems that, there are positions where both the interac- 
K Orce and moment acting on the leading ship tend to 
ам the leading ship towards the tracking ship. Such posi- 
tions seem to be at the longitudinal οι ancestro 524 τε, 
ο. τὸ and 300 ft. 

An interesting effect was observed when the 
ШОО Л Таз ха! distance is +200 ft and the lateral distance 
EE, Figure 44. Although the interaction force-Y and 
moment-N had both positive values they did not affect the 
l ra] distance (B) throughout the time of the run (2 
minutes roughly) because the interaction MCCC eeunter- 
ШЕсегсі each other, due EMEC hace hat Ehe positive inter- 
action moment was producing a yaw positive angle such as 
и соцпсекасЕ Еће positive interaction force. The the Ya 
velocity of the leading ship was zero in this case and so 
the lateral distance could be kept constant. At this point 
for the ease of the manual control the input voltage into 
meee steering control block and into the propeller speed 
blcok is changed as Figure 46 shows for ship model No l and 


Figure 47 shows for ship model No 2. 
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Figure 47. Rudder and propeller input commands diagrams 

for ship model No 2. 

Those analog configurations shown on Figures 
46 and 47 for both ship models respectively replace the up 
meme now used configuration. This is done primarily because 
One can have negative or positive rudder and propeller 
speed changes without necessarily changing the switch posi- 
mom. This Obviously gives operative flexibility for the 
operator making the control of the ship more realistic 


besides the necessary required sensitivity of the rudder 
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meee! and propeller. Then for zero rudder or propeller 
Speed change either the corresponding switch has to be open 
or the corresponding potentiometer has to be set at 0.5000 
Value. For example for the rudder angle will be 5 turns 
ID to 20° and 5 turns 0 to —20° which оссо» hom. 

Ж лКеулзе 1Е will be 5 turns for 0 to + 30 RPM 
hence it will correspond to 6 RPM per turn. For the ease 

С рРегасіоп the potentiometer values for P437, Р406, Р416, 
ШОО should be set in "RESET" mode having in mind that 100.0 
ШЕР соггесропав to 20° rudder and 30 RPM of propeller 
Speed respectively. This avoids the confusion of converting 
in "POTSET" mode the potentiometer values to be set in 
υπο volts since now the 0.5000 position of the potentio- 
meters corresponds to zero volts. 


b. Runs with different propeller speed between 
cn ME ο во and Ship B(tracking ship) 


Hwo Kinds Of runs were made, one with the track- 
πι» have a speed of +10 RPM greater than the speed of 
EN nots of the leading ship. In the second run the differ- 
ence in propeller speed was 5 RPM. Note that no controls 
were applied. The following assumption was made: both ships 
Sabe placed at any desired initial position at the begin- 
ning of the run. This was done because no complete run with 
the tracking ship passing the leading ship was possible to 
be made primarily for two reasons: 


(1) Saturation effects of the analog computer 
amplifiers. 


Штета сто effects can be applied only 
- for the lateral distance being between 
ӨШ ЕЕ апа 00 ТЕ. 
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This assumption holds for Phase II also. A matter of fu- 
ture investigation would be a rescaling of the analog com- 
puter unit to avoid the saturation effects. Also more 

а should be available for the interaction effects to 

be applied for lateral distances outside, the range of 50 
to 160 ft. Nevertheless the tracking ship was placed at 
different initial longitudinal distances from —524 ft to 
fee tt. All the responses of the leading ship for Phase 

I were calculated in terms of the perturbated parameters 

of the leading ship PY, N, Ý t ËB V S. From the observa- 
tions of the obtained responses, Figure 48 to Figure 64, 

it can be seen that during the approach and the departure 
of the tracking ship the leading ship tends to reduce the 
Шигега! distance, B. So a possibility of collision appears. 
More specifically the tracking ship starts at —524 ft longi- 
tudinal distance with 10 RPM propeller speed greater than 
that of 15 knots and tries to overtake the leading ship, 
Figure 43. Due to interaction moment N the πο ship 
yaws to a positive angle and so tends to reduce the lateral 
ЕЕ Еапсе B, although the interaction force Y is negative. 
τε 43, the next run, the tracking ship starts to 

aa reakezche leading ship at a longitudinal distance of 
—400 ft with 10 RPM propeller speed greater than that of 
mae leading ship. It can be seen that originally the yaw 
angle is positive due to the positive interaction moment 


and the lateral distance is decreasing again. This in fact 
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is happening because the bow of the leading ship yaws to— 
wards the oncoming tracking ship. When the longitudinal 
distance is reduced to —200 ft the leading ship starts to 
yaw at a negative angle; this is because both the Y-force 
апа N-moment changed sign, hence the lateral distance De 一 
meomunchnanged for about 15 sec at 50 ft approximately, 
Starts to increase toward 100 ft when the longitudinal 
dastance is approximately 160 ft. 

In the next shown run, Figure 50, similar analy- 
Sis can be made except from the fact that the tracking ship 
ШЕШ starts at initial position at е —300 ft and Vs 
ЕЕ ГІ. The lateral distance being 50 ft for about 15 sec 
Starts to increase when the longitudinal distance is approxi- 
Шету —200 ft. 

Figure 51 shows the response Orewa.) tara |. 
posstion of the leading ship at (0,0) c Ll eni r ackrking ship 
at (0,50 ft). The tracking ship possesses a 15 RMP greater 
propeller speed than that of 15 knots. The ig ship 
immediately starts to yaw negatively and thus increases the 
lateral distance drastically. 

Simebtanites 52, wiere iunttialily a = +160 ft and 
В = 70 ft it is seen that the lateral distance increases 
а | the longitudinal distance is +200 ft. Then due to 
ehe interaetion effects the leading ship reduces the 
lateral distance, since a positive Ve Velocity Пав been 


created. 
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Figures 53, 54, 55 show the response of the 
leading ship when a and ß were initially set at +300 ft, 
feo rt, +524 ft, and 70 ft respectively, at 5 RPM differ— 
ence in propeller speed. It mee nar the best position 
for two ships of the same size is exactly abeam when along- 
side. The increase of the lateral distance when the longi- 
tudinal distance is between approximately —200 ft and +200 
ft implies that the leading ship yaws. Hence control has 
ре applied for it to keep a constant course. It is also 
ΜΙ του that the approach and departure of the tracking 
ship implies the use of control on the leading ship in order 
τς tO avoid collision. 

Figures 56 till 64 Show the linear response of 
ene leading ship for different initial positions of the 
τής Ship with a difference of 5 RPM in propeller speed. 
Pe weracking Ship overtakes the leading ship, which has a 
15 knots speed. A similar analysis can be made for these 
responses as was done for the previously mentioned runs. 

ШЕ Obtained Responses for Phase II 

Computer Program VI shows the source deck used for 
Phase II. In Phase II interaction forces and moments are 
σα On beth ships. The circuitry of both Figures 33a 
and 33b is used. It should be mentioned here that the 
digital machine computes ШЕШЕСІ Ship the inter- 
action force and moment corresponding to those longitudinal 
distance, a, and lateral distance, 8, which the analog com— 


puter passes through the trunk lines to the digital computer. 
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Eug the same cycle of computation the digital computer 
inverts the sign of longitudinal distance, a, and for the 
same lateral distance, B, gives another pair of Y-force 
and N-moment to the analog computer, which are the inter- 
Exon effects applied to the tracking ship. 


a. Stationary Runs (i.e. both ships have same 
propeller speed) | 


Figures 65 to 75 show the linear response of 
both the leading ship and.the tracking ship alongside at 
the same speed of 15 knots at various longitudinal and 
lateral distances respectively. The responses were calcu- 
lated in terms of the perturbated parameters Va Y Nr Y pr 


В ὑ , Un where the indices A and B mean the leading 


ο 
ship A ama tracking Ship (ship B) respectively. 
Depending upon the relative position of the ships there 
are cases of attraction between them and repulsion between 
them. An interesting case was observed at the exactly 
abeam position, while the two ships were alongside ato 70 
ft lateral separation distance. The lateral distance was 
Rept constant at 70 ΟΠΠ ΕΠΕ time of recording 


although both ships were yawing at negative angles. 


PrzzRuns with Different Propeller Speed between 
ο ο Tesding Ship) and Ship B (Tracking Ship) 


The assumption made for Phase I that both ships 
can possess any desired initial positions with respect to 
the (ху) NC СО ОУ Тс ОЕ рһазе ІІ, since 
again the tracking ship was never able to completely over— 


make the leading ship during the progress of a single run. 


“4. 
«4 
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The obtained responses were calculated in terms of the 


perturbated parameters Var Y N 


Toe 7 A, o F Un 


where the indices A and B mean ship A and ship B respec- 


AS 


tively. With no controls applied on both ships the 
tracking ship was chosen to have in the first case 10 RPM 
propeller speed and in the second case 5 RPM propeller 
speed greater than the speed of 15 knots, which was the 
leading ship's speed. 

Iii? 76 tne initial position of the track- 
ing ship is at (-524 ft, 70 ft) and that of the leading 
AE (070). Obviously the longitudinal distance de- 
creases since the tracking ship has a greater speed. The 
lateral distance, decreases also since the Yo, velocity orf 
the leading ship is positive and greater than the Yo, veloc— 
о гле tracking ship, which is also positive. This 
means that the leading ship is pulled towards the tracking 
Ship. It is seen also that the leading ship yaws to posi- 
tive angles as well as the tracking ship but since the 
leading ship's yaw angles are greater than those of the 
tracking ship the rn of the leading ship comes towards 
БЫ ром об the tracking ship. 

n ο P у шоља Ју the longitudinal 
distance, a, and the lateral distance, ß, were -400 ft and 
HH = срескимеју the lateral distance, B, was decreas- 
L селе long+tudinal distance became roughly —250 Ees 


'ὗ τι σ1πΕ both the interaction force and moment on the 
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ESsogng shrp changed Sign. This can be observed in Figure 
78 where initially a and ß were —300 ft and 50 ft respec- 
tively. The above mentioned change of Sign results in a 
Megative leading ship's yaw angle, while the tracking ship 
continues to have positive yaw angle. Tho lateral distance 
ο. ος το increase from the instant at which the longitudi- 
Waai Stance is roughly —250 ft. 

Iuris 79 о and B were КОО О πι 
and 50 ft respectively. It can be seen that the leading 
ship yaws to negative yaw angles while the tracking ship 
maintains zero yaw angle. Although this is happening the 
I velocity is negative and the Yo, almost zero. Hence 
the leading ship is pushed away from the tracking ship and 
the lateral distance is increased. 

The same analysis can be done for the response 
Seer igure 50a where a and В were initially —160 ft ana 70 
ft respectively as well as for the response of- Figure 80b 
wijere initially was set a = ~160 ft and 8 = 50 ft. 

πο τας 8l shows the response obtained for the 
exactly abcam position of the two T being alongside. 
It is seen although both ships are yawing excessively to 
megative angles the lateral distance is slightly increas— 
ing. Hence this is the best position for station keeping. 

Figure 82 shows the obtained response for a 
and ß being initially +160 ft and 100 ft respectively. 


EM Seon that the lateral distance B starts to decrease 
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when the longitudinal distance becomes roughly +200 £t, 
while the leading ship yaws slightly negatively and the 
tracking ship yaws excessively negatively. These facts 
bring the bow of the leading ship towards the tracking 
ship's stern. At +250 ft longitudinal separation between 
midships the leading ship reverses direction of yawing 
While the tracking ship continues to yaw negatively. 

This can be seen more easily in the response 
 Елапуе 83 where initially о and B set at +200 ft and 
100 ft respectively. Due to these opposite yaw angle 
directions the sterns of both ships are pulled toward each 
Other while the lateral distance is decreasing. In Figure 
Sec and 8 initially were set to be equal +300 ft and 70 
ft respectively. The lateral distance is decreasing due 
ehe Tact that the leading ship's Yo, ο πσ ο στο posi- 
tive and greater than that of the tracking Sipe А150 
the leading ship yaws positively while the tracking ship 
ШЕМЕ Slightly negatively. This brings leading ship's 
stern towards the tracking ship. Note that at roughly 
1350 ft longitudinal separation the tracking ship reverses 
Me direction of yawing. This tends to bring the bow of 
the tracking ship towards the stern of the leading ship. 

Figures 85 and 86 show the responses for 
ss jo (Trung a and B 4400 ft and +524 ft, 70 ft res- 
Peetively. The lateral distance, 8, is increasing although 


both ships yaw positively. 
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Figures 87 to 97 show the responses obtained 
K Піттетспіе initial positions of the ships. The analy- 
sis of these curves is analogous to the previous one, 
me difference being that the tracking ship has now a +5 
en eater propeller speed than that of 15 knots. Note 
that for the smaller difference in propeller speed there 
is a smaller rate of change of the lateral separation 
distance although both ships yaw to greater:‘values of yaw 
angles since they have adequate time for that. 

a Manta ll Control is Applied in Phase II 

From the previous obtained responses it 15 
obvious that control must be applied in order to keep a 
desired course as well as a lateral separation distance, 
То τι now for lateral distance greater than 100 feet 
or less than 50 feet the digital programming was made such 
that for these cases a zero force-Y and moment-N was given 
pomehecanalog computer. This fact made impossible for the 
tracking ship to completely overtake the ns ship dur- 
ШОШО one run. At this point the following approxima- 
tion was made to override this difficulty, namely the 
digital programming, Computer Program VII, was made such 
that the digital computer gives (for values of ß outside 
the table limits) the values of Y-force and N-moment cor- 
responding to the extreme values of ß, i.e. 50 ft and 100 
ft. Obviously ME ΠΠ 5-Ε the interaction effects 
would be more pronounced for ß less than 50 ft and less 


pronounced for ß greater than 100 ft. But witn this 
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approximation it was possible to get the response of Fig- 


ure 98 with control manually applied in terms of the рек— 


turbated parameters Wns АКД, Nar Yo,’ б, р, ARA, Us where 
End cB indices mean ship A and ship B respectively. The 
Pragekıng ship originally placed at = A τὅια Ver ο πα 


Nas a 10 RPM greater propeller speed than the leading ship. 
These responses do not give the actual needed 
rudder to maintain station and course unless the operator 
1S practiced enough to know the behavior of both ships and 
hence is familiar with the amount of rudder actually needed. 
Шемен ле1ле55, these responses support the already known 
fact (from experience at sea) that underway replenishment 
1S a controllable operation although experienced helmsmen 


are needed on both ships of the UNREP. 


IPS ONCEUSIONS — REMARKS — SUMMARY 

Before any conclusion is derived it must be mentioned 
here that this study was done for the underway replenish- 
ment at sea of two equal sized merchant ships, namely 
MARINER class ships. The study also assumes calm water 
conditions which consequently assumes motion of ships in 
3-degrees of freedom only. The hybrid simulation was done 
Pony for the no control plant although some remarks 
can be deduced for the nature of rudder (control surface) 
needed. Most of the observed results of this simulation 
study are already known from experience of real replenish- 


ment at sea operations. But this agreement assures that 
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this simulation was done correctly and gives the neces— 
es background for the investigation and solution of 
the control problem. 

During the approach and departure of the tracking 
Ship there is a risk of collision Since at these posi- 
tions the lateral separation distance is decreased due 
to the interaction effects created from the pressure dis- 
pm Lion fields of the hulls of both ships. 

To be more specific, with the leading ship placed at 
(0,0) and having a 15 knots speed, i.e. the same speed 
with that of the space coordinate system, also with the 
metering Ship placed at (—524 ft, B) and having a An 
Perimogreater propeller speed than that of the 15 knots 
Ship the following can be deduced. 

a. During the approach the stern of the leading 
ship and the bow of the tracking ship tend 
kesattract, while the lateral distance is 
decreased until the longitudinal separation 
distance" το лей ТЕ, Also both ships 
yaw in the same direction, but at different 
rates and values depending on their speed 
difference. 

b. For a decrease in longitudinal separation 
ПЕЕШЕ 250 ft πο Ζει O l e lateral зерака— 
tion increases although the ships yaw to 
different directions and values. 

ЕК2 его шова оства нето апсе, 1.е. exactly 
abeam position and 50 ft lateral distance 
both ships yaw negatively by the same amount 
while the lateral separation distance slightly 
increases throughout the time of recording. 

d. With both ships having the same speed of 15 
knots and being exactly abeam there is little 
change of the lateral distance throughout 


the time of recording although both yaw 
negatively. 


195 


e. During the qeparture from the exactly abeam 
position to until the longitudinal distance 
ο ο ау 0. EE. 


ООШ σου the lateral distance takes 
puce 


КО nt trolly eee bow OL the Leading ship 
tenäs to be attracted to the stern of 
the tracking ship, while both ships yaw 
negatively at different angles. 

DE toroughb t250 ft longitudinal separa- 
tion distance the leading ship changes 
direction of yawing while the tracking 
ship Maintains negative yawing. This 
επ τ: ehe storn of the tracking 
ship towards the leading ship. 

f. For longitudinal separation distances greater 

"η τοπ, "τετ, che lateral distance 

increases while both ships yaw positively 

ΞΕ θε! πα rates. 

One basic assumption was carried out throughout the 
UNREP simulation. That is both ships could be initially 
placed at any desired position. 

From the above discussion it is seen that the best 
Pesmeron £Or Station keeping and collision avoiding for 
ps tot ле сале Size 1S while keeping the same speed 
exactly abeam, although a certain amount of rudder has to 
be applied on both ships for course keeping. This is in 
agreement with reference [1]. It is obvious that during 
the departure and approach the rudder to be applied must 
pounteract both the interaction moment as well as create 
Иа естене το counteract both the attraction 


force and the rudder force. By these means collision 


should be avoided provided that there is enough initial 


194 





lateral separation distance so the time elapsed between 
the instant action iS applied by the operator (helmsman) 
and rudder action takes place is adequate. It can be 
seen also that during the approach and departure both 
ships yaw from positive to negative angies and vice versa. 
This implies that іп а short period of time the operator 
has to apply negative to positive rudder angle provided 
EDUstbe operator has been trained enough so that he knows 
the exact time at which rudder has to be applied. Ob 一 
viously this is not easy. 

Thrs study probably would be more effective using the 
interactive graphics available on modern computer systems. 

A matter of future investigation Should be the imple-- 
mentation of automatic controls in the hybrid simulation 


il as the inclusion of sea states. 
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πι pu TD APO OPTOATIOSSOP THE -CONTROL PROBLEM 


The leading ship usually is responsible for course 
keeping and the tracking ship is responsible for both 
course and station keeping. These requirements imply two 
ШОШ TOCIL controls for the course and the station or dis- 


tance keeping. 





Steering Er 
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Figure 99. Course keeping loop. 


COURSE KEEPING LOOP 

Figure 99 shows the control loop for course keeping. 
The steering control blocx E can be represented as a 
це 6 окаек јад function. In the section B2 the time lag 


was estimated to be of the order of 8 sec. Hence the 


Fenster funetion Should be as follows: 


where Ne the steering control gain. 


τος 


Bheccourse keeping loop may not contain the yaw angle 


reference, since during the UNREP operation course 


ЕЕ 


changing is not desired but course keeping. Hence for 
the beginning at least of the simulation with automatic 
control included the reference yaw angle Ver should be 
Бес еппа! (о zero. The ships dynamics block contains the 


mieamster function G, 2 ON which can ђе found by Laplace 


transforming equations (V—1) and solving for the above 
memiemoned ratio. It should be noted here that both the 


transfer function G, and G, have been built for the analog 


V 2 


simulation of Phase I and II. Obviously in question must 


be the order of magnitude of the feedback gains Ky and K. 


The values of Ky and K can be designed by normal feedback 


control techniques. 


DT STANCE OR STATION KEEPING LOOP 


Figure 100 shows the distance keeping loop where f... 


ene Ordered lateral separation distance and B the actual 


lateral distance between the two ships. Note that this 


оК [Smp fw 
μμ] Dy NAMES <i В 
> 








тоа 100. Distance keeping loop. 
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u rO] loop contains the ships dynamics G, including 


2 

hybrid simulation for obtaining the lateral distance. 
In question must be whether the fecdback function 

Should contain the rate of change as well as the change 


of the lateral distance. The steering control transfer 


上 Me も HO is again: 


Note that for the distance keeping loop the ships 
dynamics transfer function G, involves hybrid simulation 
since transformation of coordinates is πο ο τος the 


quantities being measured with respect to the space co- 


ordinates system. 
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